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The  compliance  calibration  for  two-dimensional  cracks  is  seldom 
attempted.  This  is  due  to  the  fact  that  a fairly  wide  selection  of 
methods  [e.q.,  finite  elements,  collocation,  complex  variable.  Inte- 
gral equations,  etc.)  are  available  to  obtain  stress  intensity  factors 
to  a good  degree  of  accuracy  and  also  due  to  difficulty  in  performing 
experiments  with  sufficient  accuracy.  For  the  three-dimensional 
problems,  in  general,  the  difficulties  are  compounded,  analytically  as 
well  at  experimentally.  However,  for  an  internally  pressurized  thick- 
wall  cylinder  with  a symmetric  internal  surface  flaw  located  in  the 
axial  direction,  the  evaluation  of  compliance  can  be  accomplished. 

This  has  been  done  by  proper  definition  of  compliance  and  data  reduction. 

In  fact,  a fatigue  crack  originating  at  the  bore  surface  tends  to  be 
orientated  with  its  length  In  the  axial  direction  and  depth  in  the 
radial  direction. 

The  theoretical  basis  for  the  compliance  test  using  the  linear 
theory  of  elasticity  is  given  in  Section  one. 

With  a starter  notch  of  a semi-circular  cross-section  an  extensive 
set  of  tests  were  performed  as  explained  In  Section  two.  In  the  third 
section,  the  method  of  numerical  data  analysis  is  outlined,  and  the 
stress  intensity  factors  are  computed.  In  Section  four,  we  simulate  the 
actual  crack  by  a simpler  geometry  and  using  the  method  of  three-dimen- 
sional collocation,  the  stress  Intensity  factors  are  computed.  The  exper- 
imental results  compare  well  with  the  theoretical  results  obtained.  To  the 
authors'  knowledge,  such  a calibration  does  not  exist  in  the  literature. 
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SECT  I CM  i - ANALYSIS  OF  COMPLIANCE  CAI  IBRAT1  (IN 

The  volumetric  compliance  method  for  a thin,  straight-fronted 
radial  notch  has  Seen  discussed  in  [1],  For  the  present  problem  of  a 
symmetric  surface  flaw  located  on  an  axial  plane  shown  In  Flqure  1,  it 
is  ciear  that  strain  energy  release  rates  and  stress  Intensity  factors 
are  functions  of  the  axial  coordinate  x^.  However,  If  we  restrict  our 
attention  to  the  plane  x^O  only,  the  concepts  of  volumetric  compliance 
can  be  extended  to  compute  the  stress  intensity  factors  at  the  deepest 
point  of  the  surface  flaw,  i.e.  at  x^=0.  Let  the  following  quantities 
be  defines  al  AV.  is  the  change  In  the  Internal  volume  per  unit 
length  In  sue  to  the  internal  pressure  p;  G Is  the  crack  extension 
force  (strain  energy  release  rate);  C^.  = 1/p  Is  the  volumetric 

compliance,  and  b is  the  crack  depth.  Then, 

, o dC 

G ' - ' at  xn“°  ' 0 ) 

2 db  3 


The  stress  Intensity  factor  at  Xj=0  , denoted  by  Is  , 

2 

K 2 = — T Wf/P)  • at  x,*0  . (2) 

1 (1-v2)  2(1-/)  db  1 3 

The  above  relation  is  a consequence  of  a 'plane  strain1  condition  at 

x3=0,  [2].  As  Indicated  In  (2),  the  compliance  calibration  can  be 

'UNDERWOOD,  J.H.,  LASSELLE,  R.R.,  SCANLON,  R.D.  and  HUSSAIN,  M.A., 

"A  Compliance  K Calibration  For  A Pressurized  Thick-Wall  Cylinder 
With  A Radial  Crack"  Engineering  Fracture  Mechanics,  Vol.  4,  1976, 
pp.  231-244. 

2EMERY,  A.F.  and  SMITH,  F.W. , "Some  Basic  Properties  of  Penny-Shaped 
Cracks"  Mathematlka,  Vol.  13,  1966,  pp.  172-180. 


accomplished  at  a given  value  of  pressure  by  determining  the  change  in 

the  internal  volume  of  the  cylinder  (at  X3=0)  for  a series  of  surface 

Haws  of  different  depths  and  performing  the  indicated  differentiation 

However,  the  measurement  which  can  easily  be  made  is  the  change  in  the 

outside  volume  cf  the  cylinder  per  unit  length  at  x^O  (equivalent 

to  the  change  in  cross-sectional  area  at  x^O).  In  the  sequel  we  will 

prove  that  the  derivative  of  internal  and  external  volume  changes  for 

a given  pressure  is  identical.  AVQ  and  AVm  are  the  changes  in  the  out 

side  volume  and  material  volume  per  unit  length  at  x^O.  Hence 

AVq  = AV.  + AVm  (3) 

If  <5V  is  the  total  change  in  the  material  volume  V bounded  by 

cylindrical  surfaces  and  between  cross-sections  S^.:  and 

S*  x and  exterior  to  the  upper  and  lower  surfaces  of  the  surface 
63 

flaw  and  , Figure  1,  we  have 

6V  = /v  e dV  , (4) 


where  e - ls  dilation.  The  generalized  Hooke's 
law  *or  a homogeneous  isotropic  body  can  be  written  (using  the  usual 
notations)  in  the  following  form:  (see  e.g.  [3]) 


Tij = uije  + 


(ij  = 1,  2,  3) 


(5) 


Contraction  gives: 


Hence 


x..  = (3A  + 2y)e  . 

* j 


5V 


1 

3,\+2y 


x . . 
1 1 


d V 


(6) 

(7) 


“'30K0LNIK0FF,  1.5. , "Mathematical  Theory  of  Elasticity"  McGraw  Hill 
3ook  Company,  Mew  York,  1956,  p.  66. 
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Using  the  equation  of  equilibrium 

= 0 , (i ,k  = 1 , 2,  3)  , 

one  can  show: 


fvVWv(Tik  V’*  dV 


Tne  divergence  theorem  yields 

J (t.il  x-)>t  dV  = J v t x.  dS  = / T.,  x.  dS 
Jy  Ik  o x 's  k lk  3 S 1 3 


(8) 

(9) 


(10) 


where  S 's  a surface  forming  the  complete  boundary  of  the  volume  V 

v 

and  o,-.'  the  direction  cosines  of  the  outer  normal  v of  S and 

v 

are  .no  ccrponents  of  the  traction  T acting  on  S.  Hence  we  have 
from  (9)  jr.o  ;10) 

jjt|jd,.|sT,  xjds  on 

and  contraction  gives 

Jv  'ii  M ' /s  T1  x,  dS  = /s  T . rds  (12) 


The  surface  5 in  Figure  1 is  formed  by  six  surfaces  denoted  by  , 
So , • • • , S6  then 


,v  ,iv 

111  2 


l ' 2 i 3A+2ii  nil  (T>)h  • r 


dS 


(13) 


k+2u  n»)  ~n 

0 

where  (T)n  are  tractions  on  the  surface  Sn>  From  the  symmetry  we 

v v v v 

have  (T)^  = (T)^  and  (T)g  s (T)g  and  the  outer  normals  are  in  opposite 
directions  for  each  pair  (Sj,  S^)  and  (Sg,  Sg).  The  last  four 
integrals, having  possibly  the  crack  depth  b as  a parameter,  vanish. 


This  completes  the  proof  that  AVm  is  independent  of  b. 
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CIGURE  2.  Schematic  diagram  showing  the  specimen  for  compliance  test 


(U) 


For  O'* 
where  c 

and  free- 


profiurued  cylinder  we  have  (T) -j  = -er  p , (T)2  = 0, 
,;'  i . vector  ir,  the  radial  direction.  Thus 


dij  db 


or. 


The  stress  intensity  factor  then  can  pe  written  in  terms  of  the  out- 
ci:‘e  vr-'.itc  change  of  the  cylinder  per  unit  length 


■ !>  - 


!2(l-vW 


d(AV0/p) 

— ar~ 


1/2 


at  xyO 


(16) 


Si  ■ ;:•'!#  i .invU  Ti’Si' 

The  spec’r'en,  made  from  a steel  cannon  tube,  was  a cylinder  with 
a 7.1  Inch  smooth  bore  and  a 14.25  inch  outside  diameter,  giving  a 
diameter  ratio  of  2.007  and  a wall  thickness  of  3.575  Inches.  The 
length  was  30  inches  with  a recess  1.5  Inch  long  and  7.3  inch  diameter 
at  each  end  for  ‘ h<  pressure  seal. 

if  was  tested  as  a vertical,  hollow  cyl inder,  internally  pressurized 
on  a supporting  mandrel  so  that  the  end  loads  were  supported  by  the 
internal  mandrel,  giving  an  open  end  condition.  (See  Figure  2.)  The  cyl- 
inder ends  were  sealed  with  threaded  closures  on  the  mandrel  so  that 
there  were  no  axial  stresses  produced  in  the  specimen  due  to  pressure 
end  loads.  The  pressure  fluid  was  a synthetic  Instrument  oil  which 
remains  fluid  at  pressures  In  excess  of  those  used.  The  pressure  was 
fully  exerted  in  the  initial  notch  cavity  when  the  cylinder  was  pres- 
surized. 
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Table  1.  Circumferential  Strains  (in. /in. 

Pressure  for  the  Successive  Crac 
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The  wall  was  initially  notched  with  a semi-circle  of  1/4  inch 
radius  cut  at  the  bore  by  electrical  discharge  machining  (EDM).  The 
notch  was  1/2  inch  long  in  the  axial  direction,  1/4  inch  deep  in  the 
radial  direction  and  0.030  inch  wide  In  the  circumferential  direction 
at  the  midsection  of  the  specimen.  A fatigue  crack  was  grown  from  this 
starter  notch  by  hydraulic  pressure  cycling  from  zero  to  48,000  psl 
until  it  reached  each  successive  desired  depth,  as  determined  by  ultra- 
sonic crack-depth  measurement.  The  circumferential  strains  were  then 
read  on  the  outside  circumference  with  bonded  foil  strain  gages  at  14 
angular  positions  relative  to  the  notch  location. 

The  strain  gages  were  1/4  inch  gage  length  and  were  oriented  In 
tne  circumferential  direction  in  a ring  at  specimen  midlength.  Two 
gages  were  located  on  a line  directly  over  the  notch,  one  at  1 inch 
above  and  one  at  1 Inch  below  the  cylinder  midlength,  so  as  to  permit 
measurement  of  crack  depth  at  its  deepest  point  by  placement  of  an 
ultrasonic  normal  probe  directly  over  the  notch.  Other  gages  were 
located  around  the  midlength  circumference  at  angles  from  5°  to  180° 
on  either  side  of  the  notcn  line,  as  listed  in  Table  1. 

When  measuring  the  strains  at  each  crack  depth, the  pressure  was 
increased  to  16000,  32000  and  48000  psl  and  then  decreased  to  40000, 
32000,  24000,  16000,  8000,  4000  and  zero.  The  pressures  were  controlled 
within  + U.  The  readings  from  all  gages  were  taken  in  sequence,  with 
the  cylinder  held  at  pressure,  using  a 24  channel  switching  unit  and  an 
SR-4  strain  Indicator  located  outside  the  test  cell.  An  identical 
strain  gage, bonded  to  an  unstrained  piece  of  the  same  steel  lying  on  the 


9 


test  specimen  provi  ’ed  temperature  compensation  during  the  experiments. 
If  the  final  strain  reading  at  zero  pressure  differed  from  the  Initial 
zero  pressure  reading  by  more  than  20  microinches  per  inch  at  any  crack 
depth  the  readings  were  repeated  to  eliminate  measurement  errors  as  much 
as  possible. 

Strain  gage  data  for  the  48000  psi  pressure  level  are  listed  for 
the  successive  crack  depths  in  Table  1.  The  gage  positions  are  the 
angular  locations  on  the  0.0.  from  a point  directly  over  the  initial 
notch. 

figure  3 shows  the  final  fracture  surface  of  this  fatigue  crack. 

The  crack  profile  at  successive  crack  depths  is  reproduced  in 
Figure  4.  Tne  corresponding  number  of  cycles  to  reach  each  depth  is 
given  in  Table  1.  Figures  3 and  4 show  that  the  crack,  propagated  from 
the  initial  semi-circular  notch,  remained  a semi-circle  until  nearly 
2-inch  depth  was  reached.  After  that  it  became  part-circular  or  part- 
elliptical.  It  finally  acquired  the  shape  of  a partially  imbedded 
ellipse,  when  the  intersection  of  the  crack  front  with  the  bore  surface 
became  nearly  stationary  while  the  ellipse  deepened  and  lengthened.  In 
these  latter  stages  of  growth,  after  a 2-in.  depth,  the  rate  of  propaga- 
tion was  quite  fast.  Unfortunately  we  did  not  obtain  as  complete  strain 
data  as  would  be  required  for  an  accurate  compliance  analysis  for  the 
greater  depths.  The  data  is  therefore  analyzed  for  the  semi-circular 
shape  up  to  the  2-1nch  depth  only.  In  Figure  3 one  can  see  the  evidence 
of  a “free  surface  shear-lip"  formation  which  fans  out  at  about  20° 
from  the  inner  wall  on  either  side  of  the  initial  notch. 
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FIGURE  3.  The  final  fracture  surface  of  a Fatigue  Test  Specimen 


FIGURE  4.  Crack  profiles  at  successive  crack  depths. 


SECTION  3 - DATA  /WJALYSIS 

by  considering  only  the  cross-section  of  the  cylinder  through  the 
strain  gages  (i.e.,  X3=0),the  numerical  problem  is  stated  as  follows. 
Given  the  tangential  strain  at  selected  points  on  the  circumference  of 
a circle, we  wish  to  determine  numerically  the  change  in  cross-sectional 
area  (equivalently  AVq  , the  change  in  external  volume  per  unit 
length  at  x^-0) . 

This  can  be  accomplished  in  two  ways.  First, we  can  integrate  the 
tangential  strain  to  obtain  the  change  in  the  length  of  the  perimeter 
AP  = d<i>  , at  Xj=0  , (17) 

and  assuming  the  distorted  figure  remains  circular,  we  have  the  change 
in  cross-section  area  AAq'. 

M0  = 1 (2PAP  + AP2)  . (18) 

A more  rigorous  approach  is  to  compute 

AAn  = J2"  fr2+ur(r><)>)  ,.drd4,  (19) 

0 0 r2 

From  strain  displacement  relations, we  have  (using  cylindrical 
coordinates): 

3u  (r  ,4,) 

ur(r2,<(i)  = r2ed,(r2  ’^ ) ^ ’ (20) 

where  u , u are  radial  and  tangential  displacements.  It  can  be 
r 

shown  by  eigenfunction  expansions  that  the  second  term  at  the  right  of 
(20),  for  the  present  problem,  is  even  in  $ and  periodic  in  2ir  . 
Neglecting  second  order  terms  (linear  elasticity),  It  is  seen  that  this 
term  makes  no  contribution  to  the  Integral  in  (19).  Hence, essentially 
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we  have  (r2*<*>)  5 r2e4>^r2 ^or  the  inte9rat''on  of  (19)-  It  was  found 
that  the  evaluatlonsof  Mg, using  the  two  methods, agree  with  each  other 
to  four  significant  figures. 

For  the  data  analysis, our  choice  of  an  angular  function  was  an 


interpolating  periodic  cubic  spline  [4], 

Definition  of  cubic  spline: 

Given  an  interval  , a£X£G  , (21) 

a mesh  on  the  interval, 

A ; a = x0  < x1  < *•••  < xN  = B , (22) 

and  associated  set  of  ordinates  (data  points)  , 

Y : yQ,  y1 , •••  , yN  (23) 

then  the  cubic  spline  satisfies: 

$^(Y;  x)  € C2  on  [a,  £l]  > 

Sfi(Y,  Xj)  = yj  (j=0,  1.  •••  , N)  , (24) 

and  is  coincident  with  a cubic  on  each  interval, 

*j-l  i x < Xj  (j»l,  2.  •••  , N)  . (25) 

If, in  addition, the  derivatives 

SA<P)  (a+)  “  4  Sa<P)  (p=0’  11  Z)  (26) 

the  spline  is  said  to  be  periodic  with  period  ( B-a)  . 

If  alternatively 

sa(Y;  Xj)  = yj  + Ej  , (o=0,  1.  — , N)  (27) 


with  the  cj  subject  to  some  minimizing  constraint;  it  is  said  to  be 
an  approximating  rather  than  interpolating  spline. 

4AHL8ERG,  J.  H.,  NILSON,  E.  N.,  WALSH,  J.  L.,  "The  Theory  of  Splines 

and  Their  Applications,"  Academic  Press,  New  York,  1967. 
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As  described  in  the  previous  section  the  data  was  taken  at  14  angular 
position  for  six  crack  depths  and  at  nine  pressures.  Discounting  the 
data  for  crack  depths  greater  than  2 in.  (due  to  gage  and  tube  failures) 
we  have  756  data  points  for  numerical  analysis.  To  indicate  the  smooth 
nature  of  the  strain  data  and  interpolating  functions  we  show  in 
Figures  5a  - 5e  the  results  for  various  crack  depths  at  48,000  p.s.i. 

Or.ce  the  spline  function  has  been  found  the  changes  In  the  cross- 
sectional  area  are  computed  by  the  methods  indicated  above.  The  partial 
results  ere  presented  in  Table  2.  It  should  be  noted  that  the  first  row 
in  Table  2.  for  zero  crack  depth,  is  not  empirical  data  but  extrapolated 
values  and  the  set  of  data  corresponding  to  starter  notch  (b-. 25)  have 
been  eliminated.  The  stress  intensity  factors  now  have  to  be  computed 
from  (16) 

1/2 

, at  x3=0  , (16) 


K,/P  = 


E dtoyp)' 


2(l-v2) 


db 


where  AVq  Is  the  change  in  the  outside  volume  of  the  cylinder  per  unit 
length  at  x^=0  and  thus  is  numerically  equivalent  to  AAq  and  also  from 
E d(AV0/p)_'1/2 


K^/p  = 


db* 


, at  x3=0  , 


(28) 


where  b*  = b + --  (— ) Is  an  effective  crack  depth  suggested  by 
2 it  «y 

Irwin  [5],  (0  < h < 1 , = yield  stress)  to  correct  for  the  effect  of 

plastic  deformation  at  the  crack  boundary. 


5 1 RW l N , G.  R. , "Structural  Aspect  of  Brittle  Fracture1',  Applied 
Materials  Research,  Vol.  3,  65,  1964. 
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To  compote  the  derivative  of  AV^ ( b)  we  used  a non-periodic  approximating 
cubic  spline  with  an  algorithm  which  minimized 

bn  n ? 

' £3" (b) } db  + X T (S  (b.)  -AV  . y (29) 

'0  ’>  i=i  A 1 01 

where  X is  introduced  to  allow  us  to  strike  a balance  between  the 
amount  of  smoothing  desired  vs.  maintaining  the  integrity  of  data.  As 
seen  from  (28)  we  must  of  necessity  evaluate  b*  and  K^*/p  » for  a 
given  value  of  h,  by  an  iterative  process  which  is  stable.  As  a guide 
to  selection  of  \ we  consider  the  type  of  expression  common  in 
fracture  mechanics  , = (constant)  and  see  that  d^K^ /db2  is 

strictly  monotonic  increasing.  A value  of  X=1 0 proved  quite  reasonable 
leading  to  a rapid  and  stable  convergence  of  K^*/p. 

The  results  of  such  an  analysis  for  pressures  of  40,000  and  48,000 
psi  are  shown  in  Figure  6.  in  an  ideal  linear  experiment  K^/p, 
for  different  pressure,  should  fall  on  the  same  curve.  To  partially 
account  for  plastic  deformation  at  the  crack  surface  we  have  incorporated 
the  Irwin's  correction  term,  as  explained  before.  In  Figure  6,  we  show 
the  results  for  48,000  and  40,000  psi  without  the  Irwin-correction 
(h=0)  and  with  correction  (h=l).  For  clarity  we  have  deleted  the  results 
for  intermediate  pressures  and  different  values  of  h.  In  Figures  7,  8, 
and  9 we  show  the  results  for  different  pressures,  and  h.  In  all  of 
these  figures  it  is  seen  that  the  results  are  fairly  close  to  each  other. 

SECTION  4 - RESULTS  BY  COLLOCATION  METHOD 

Because  of  the  complexity  In  the  geometry  of  the  problem  it  seems 
that  an  analytic  closed-form  solution  is  not  yet  possible.  However, 
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32000  PSI 


FIGURE  7a.  Pressure  increased  to  32  ksi . 


CRACK  DEPTH  (IN) 


FIGURE  7c.  Pressure  decreased  co  40  ksi 


CRP.CK  DEPTH  (IN) 

FIGURE  7d.  Pressure  decreased  to  32  k 


CRACK  DEPTH  (IN) 


FIGURE  9a.  Pressure  increased  to 


480Q0 


FIGURE  9b.  Pressure 


uoooo 


Pressure  decreased 


FIGURE  9d.  Pressure  decreased  to  32  ksi . 


since  we  are  working  with  relatively  shallow  surface  cracks  compared  to 
the  wall  thickness  of  the  tube  (most  of  the  fatigue  life  is  confined  to 
such  crack  growth  as  seen  from  Table  1)  we  can  effectively  simulate  the 
curved  surface  of  the  tube  by  a rectangular  block  with  a surface  flaw 
in  the  form  of  a segment  of  a circle  as  shown  in  Figure  10  and  apply  the 
method  of  collocation.  It  is  evident  from  Figures  3 and  4 that  there 
is  a 'shear  lip1  region  extending  radially  to  about  20°  on  either 
side  of  the  inner  wall  of  the  tube.  Hence  the  actual  crack  was  sim- 
ulated by  a segment  of  a circle  with  a matching  depth  and  the  center 
having  an  offset  d of  .3(radius)  of  the  circle  as  shown,  i.e.(b=a+d). 

The  thickness  (TH)  of  the  block  was  taken  equal  to  the  wall  thickness 
(i.e.  TH=Rq-Ri)  and  it  was  found  that  it  is  sufficient  to  take  half  the 
width  (HW)  and  half  the  height  (HH)  equal  to  three  times  the  radius  of 
the  crack  (HW=HH=3a).  The  loading  condition  on  the  block  corresponds 
to  the  Lamd  solution  together  with  uniform  tension  due  to  the  presence 
of  pressure  inside  the  crack: 


az/p  = 1 + 


ri2  + 
T7 


R. 


R0  (Ri+d+x) 


(30) 


Ro2 


In  the  sequel,  we  give  a brief  outline  of  the  exterior  collocation 
method,  the  details  can  be  found  in  [6]. 


6HUSSAIN,  M.  A.,  HAGGERTY,  R.  G.,  PU,  S.  L.,  and  NOBLE,  B.,  "Exterior 
Collocation  For  Three  Dimensional  Surface  Flaw",  WVT-TR-75053,  Water- 
vliet  Arsenal,  Watervliet,  N.Y.,  1975, 
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The  method  of  exterior  col  location, for  a symmetric  surface  flaw  located 
on  the  plane  z=0,  requires  construction  of  a set  of  functions  which  sat- 
isfy the  mixed  boundary  conditions  on  the  plane  z=0.  These  conditions 
are:  the  normal  stress  is  zero  inside  the  crack-surface  and  normal  dis- 

placement is  zero  outside  the  crack-surface.  Symmetry  of  the  problem 
further  requires  shearing  stresses  to  vanish  on  the  entire  plane  z=0.  For 
the  case  of  a circular  flaw  this  was  accomplished  using  three-dimensional 
Boussinesq  potentials  where  the  harmonic  functions  were  represented  by 
Kobayashi  potentials  (integral  representation  of  potential  functions) 
and  superposition  of  such  a solution  with  uniform  field  and  plane  strain 
solutions  to  satisfy  the  above  boundary  conditions.  The  results  are: 

or  = l An(n(n+1  )r'^[(l-2v)Irzl2]-(I3'Zl4)-r'1[(l-2v)Ig-zI6] 
n 

-Ena'n'3^2(n-2)(2n+l)r(n+l/2)rn/[4^vr(n+l)]}  cos  ne  , 

o9  = l An{-n(n+l)r-2[(l-2v)I1-zI2]-2vI3+r-1[(l-2v)I5-zI6] 
n 

+c„a'n'3/2(n+2)(2n+l)r(n+l/2)rn/[4/?vr(n+l)]}  cos  ne  , 

oz  = A0/S72a"3/2+£  An(-I  -zl4+ena'n'3/2(2n+l)r(n+l/2)rn/ 
n 

[/2r(n+l )]}  cos  n6  , 

tr6  = l t:n  An{n(n+1  )r_2[ ( 1 -2v) I -z I2]-nr_1  [ (1 -2v>)  1 5-z I6] 
n 

+a_n“3/2n(2n+l )r(n+l/2)rn/[4/'2\ir(ii+l )])  sin  ne  , 

Trz  = An  zt-nr'1 13+i7]  cos  ne  , 
n 

t9z  = l - enAnnzr"ll3  s1n  ne  • (31) 

n 
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Where  cn=0  for  n=0,  en=l  for  n/0  , the  summation  l is  over  n=0,  1/2, 

n 

1,  3/2,  •••,“  and  An's  are  constants  to  be  determined  by  the  collocation 
method.  We  have  used  the  cylindrical  coordinate  system  (r,e,z)  with  the 
origin  located  at  the  center  of  the  cricular  crack  of  radius  'a1,  and  I's 
are  integrals  given  below: 

11  = /”  k'2  Vkr)  Jn+3/2Cka)  k1/2e'zkdk  , 

12  = - 311/3Z  , I3  = - 3l2/3z  , I4  = - 3I3/3z  (32) 

and 

= k ljn-l^kr^  Jn+3/2^ka^k  Zkdk 
I6  = - 3I5/3z  , I?  = - 3I6/3z  (33) 


After  some  basic  integration.lt  can  be  shown  from  (31)  above, 
0 , r < a , 


oz(z=0)=  / 


(34) 


4 !♦:  ""'Jl1'1"1' ; f(4  ;•  4 


r,n  (2n+1)f(r+l/2) 
r (n+l ) 


} cosnO,  r>a. 


( 0 , for  r > a , 


2Guz(z=0)=  / 


2/27it  (l-v)a“3^2(a2-rz)1/2  l An(r/a)n  cosne  , r<a  . (35) 

( n 

Hence  it  is  seen  from  (33),  (34)  that  the  set  of  functions  selected 
satisfy  the  mixed  boundary  conditions  of  the  crack  problem.  The  stress 
intensity  factor  can  be  computed  from  equation  (34)  or  (35): 
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Table  3.  Values  of  An's  for  N=28,  u=0.3  and  a=l",  d=0.3",  TH=3.575", 
HW=3",  KH=3"  with  Linear  Constraint  A„=  -VA-cos  n 
(cos-1 (d/a)} 


r. 

r~ 

An+l/2 

0 

0.17074000 

xlO1 

- 0.28808100 

xl  02 

1 

0.21008300 

xlO3 

- 0.72235200 

xl  O-3 

2 

0.13196900 

xlO4 

- 0.21668300 

xlO4 

0.21888900 

xlO4 

- 0.15776100 

xlO4 

4 

0,78043600 

xlO3 

- 0.22554700 

xl  03 

‘a 

0.51905600 

xlO1 

0.24505900 

xlO2 

. 

4 

- C. 89767900 

xlO1 

0.17788400 

xlO1 

7 

- 0.18092500 

x 10^ 

0.16957900 

xlO1 

8 

- 0.66020300 

o 

o 

X 

- 0.28821800 

xlO'1 

9 

0.13661300 

xlO0 

- 0.43773500 

xlO"1 

10 

- 0.13656300 

xlO"1 

0.18752300 

xlO"1 

11 

- 0.89871200 

xlO'2 

0.26154100 

xlO'2 

12 

- 0.48300500 

xlO’3 

0.50154900 

xlO"4 

LI_ 

- 0.14393100 

xl0"3 

- 0.15544700 

xl  O'6 

40 


(36) 


K = K (e)  = /l  a"1  l A cosne 

11  S " 

once  the  A^'s  have  been  computed. 

In  general  it  is  not  possible  to  obtain  A^'s  in  the  infinite  series. 
An  approximate  solution,  however,  can  be  found  by  satisfying  boundary 
conditions  at  m selected  discrete  points,  called  collocation  points, 
on  the  remaining  boundaries  exterior  to  the  plane  of  the  crack.  At  each 
boundary  point  there  are  three  boundary  conditions.  From  the  3m  equa- 
tions we  can  at  most  solve  for  3m  unkown  coefficients  A 's.  To  avoid  an 

n 

illconditioned  matrix,  the  boundary  conditions  are  to  be  satisfied  in 
a least  square  sense.  Computationally,  the  least  square  procedure,  in 
obtaining  stable  results,  is  to  choose  a rather  large  value  of  m and 
solve  A^'s  from  the  N*N  system  of  'normal  equations',  which  are  derived 
by  minimizing  the  sum  of  the  squares  of  errors  at  collocation  points. 
Since  the  value  of  N is  not  known  a priori,  we  solve  the  normal  equa- 
tions for  N=1 ,2, . . .until  the  required  stability  is  achieved  in  and 
there  is  no  significant  change  in  'residuals'.  As  explained  in  [6] 
there  Is  an  additional  requirement  for  surface  flaw  problems;  namely 
the  boundary  traction  must  vanish  at  points  of  the  crack-free  surface 
intersection.  This  gives  a linear  constraint  on  An's  (see  Table  3)  and 
leads  to  zero  stress  intensity  factors  at  such  points.  With  Til,  HH,  HW, 
a,  d,  v,  NjFigure  10;  as  input  parameters  and  automatic  generation  of 
m nearly  equally  spaced  collocation  points,  a computer  program  was 

^HUSSAIN,  M.  A.,  HAGGERTY,  R.  G.,  PU,  S.  L.,  and  NOBLE,  B.,  "Exterior 
Collocation  For  Three  Dimensional  Surface  Flaw",  WVT-TR-75053,  Water- 
vl let  Arsenal,  Watervliet,  N.Y.,  1975. 
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OF  COEFF. 

of  collocation  results  for  N>15 
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■autofrettage 


written  which  calculates  An's  and  the  stress  intensity  factors  as 
briefly  described.  (The  symmetry  of  the  problem  allows  us  to  work  with 
half  of  the  geometry.) 

In  Figure  11,  we  plot  K^/p  vs  N,  up  to  N=35,  showing  the  stability 

of  the  results.  It  is  seen  that  the  results  are  stable  within  3-4%  for 

N>15.  In  Figure  12,  we  plot  K^/p  vs  6 for  N=28  with  109  collocation 

points.  It  is  seen  that  K^/p  is  maximum  at  6=0  (i.e.  at  the  deepest 

point  of  the  crack)  and  vanishes  at  6=107.46°  (i.e.  at  points  of  the 

crack-free  surface  Intersection).  In  Table  3,  we  list  A 's  for  the  above 

n 

case  showing  convergences  of  Afi  for  large  n.  Similar  computations  were 
carried  out  for  crack  depths  of  b=a+d  = .1,  .2,  .3,  .4,  .5,  .75,  1.0, 

1.3  and  2.  Inches.  The  results  are  shown  in  Figure  13.  In  the  same 
figure,  we  have  plotted  K^/p  obtained  by  compliance  results  (see 
Figure  9),  Excellent  agreement  Is  seen  between  the  experimental  and 
theoretical  results. 

CONCLUDING  REMARKS 

In  this  report  we  have  developed  a compliance  calibration  of  a 
thick-walled  circular  cylinder  with  a syjimetric  internal  surface  flaw. 
Results  obtained  agree  well,  for  small  crack  depths,  to  those  obtained 
analytically  for  a semi-circular  flaw  by  the  collocation  method. 
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APPENDIX  A 


A collocation  analysis  for  simulated  crack  geom.try  was  also 
carried  out  for  the  case  of  an  autofrettaged  tube.  Due  to  the  presence 
of  compressive  stress  at  the  inner  bore  of  the  tube,  it  was  reasonable 
to  assume  the  loading  condition  as  a superposition  of  applied  and  the 
residual  stress  fields.  Hence,  eg  (30)  was  replaced  by: 

°z/P  ' ’ + (1-aW)  { + U7dfe2  ’ 


( °y  , ,a+d+x, ...p2,  ,a2.  a2 


^[(2  log(^)+l^)-(^ 
P>3  p 


b2'A2,(a+d+x)2')(2  l03^V)/(1'b2)L 


for  (a+d+x)<p 


[ (p2/b2+p2/ (a+d+x)2 )- (2  logp/a+1 -p2/b2) 


(A-l) 


(a2/b2+a2/(a+d+x)2)/(l-a2/b2)]. 


for  (a+d+x)>p  , 

where  a=R. , b=R  , are  the  Inner  and  outer  radii,  o the  yield  stress  , 

io  y 

p is  the  pressure  and  p Is  the  Interface  depending  upon  the  percentage 
of  autofrettage  (e. g. , for  3056  autofrettage  p»a+.3(b-a).)  The  results 
of  such  a computation  for  R^3.55“  Ro=7.125"  p=48,000  psi,  oy=170,000 
psi,  and  crack  depths  b = .125”,  .25",  .5",  .75",  .9",  1",  1.25",  1.5", 
1.75",  2"  are  plotted  in  figure  A1 . It  Is  seen  that  K^/p  is  substan- 
tially lower  for  an  autofrettaged  tube  as  compared  to  a non-autofrettaged 

3K]/P 

tube  (Fig.  13).  It  Is  further  seen  that  — r — is  bounded  as 
opposed  to  a non-autofrettaged  tube  at  b=0. 

The  compliance  test  as  previously  explained  was  done  for  a 30* 
autofrettaged  tube  of  the  same  dimensions  as  before.  The  change  in  the 


47 


FIGURE  A1 . Graph  of  K /p  vs  b for  30%  autofre' 
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cross-sectional  area  was  computed  as  before  and  the  results  for  a few 
pressures  are  given  In  Table  Al.  It  is  immediately  seen  that  the  data 
is  quite  scattered  for  small  crack  depths.  This  may  be  due  to  the  fact 
that  the  notch  was  put  in  before  the  autofrettage  process  or  may  be 
due  to  relaxation  of  residual  stresses  under  cyclic  loading  in  the 
presence  of  a notch.  Neglecting  the  first  three  data  points  and  extrap- 
olating the  data  for  zero  crack  depth  to  AAQ=.357in^  the  stress  intensity 
factors  with  h=1.0  and  X=1 000  (less  smoothing)  were  computed.  The 
results  are  shown  in  Figure  Al , together  with  the  collocation  data. 

Again,  fairly  good  agreement  is  seen  for  small  crack  depths. 
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